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Possibility for the formation of stable Al atomic wire on the Si(211) surface is investigated using 
density functional theory based total energy calculations. The stable adsorption sites and the 
surface structures at various sub-monolayer coverages of Al are presented. It is found that the 
most stable and natural surface structures around one monolayer coverage is either 1 x 5 or 1 x 6 
which agrees with experimental observations. More significantly, our study revealed that unlike 
the case of Ga the formation of continuous atomic Al chains assisted by interrupted Al chains may 
be possible by controlling the experimental conditions (pressure and temperature) such that the 
Hai remains within ~ 0.3 eV from the bulk value. While the Al covered Si(211) surface may be 
metallic or semiconducting, the Ga covered Si(211) is always semiconducting in nature below 1 ML 
coverage. Also significant is the fact that Al atoms move from groove sites to the sites between 
terrace and trench atoms through lower energy channels as the coverage goes from 1/8 monolayer 
to 1/4 monolayer. 

PACS numbers: 73.20. -r, 73.21. Hb, 73.90. +f 



I. INTRODUCTION 

Technological importance of nanowire structures on sil- 
icon substrates has generated much interest in study- 
ing structural arrangements of metals on various Si 
substratesiS, In particular, metals like Al, Ga and In at 
sub-monolayer coverages on the Si(001) surface have been 
examined^ for the formation of stable atomic wires on 
surfaces. Though, stable atomic wire structures have not 
been found on the bare Si(001) surface, there has been 
some progress^Sii towards the formation of stable atomic 
wires on the patterned hydrogen terminated Si(001) sur- 
face. 

Recently, the vicinal Si surfaces have attracted much 
attention for the growth of self assembled atomic 
wires on them. The Si(211) surface 8,3 is one of 
the vicinal Si surfaces that is being widely used by 
experimentalistsiSiiiiiSii^ as a substrate for growing 
atomic wire structures. This surface can be viewed as 
a stepped arrangement of narrow (111) terraces. 

A three dimensional view of a small portion of the ideal 
Si(211) surface is shown in Fig.Q] The atoms marked T 
(called the terrace atoms) on the terrace are three-fold 
coordinated and thus have one dangling bond each; those 
on the step edge, marked E (called the edge atoms) are 
two-fold coordinated and have two dangling bonds each. 
Second layer Si atoms which are not directly bonded ei- 
ther with terrace or edge atoms are denoted as Tr (called 
the trench atoms) have one dangling bond each. The 
Si(211) surface consists of two- atom wide terraces be- 
tween terrace and edge atoms along the [111] direction. 
Two consecutive terraces are separated by steps and are 
9.4 A apart in the [111] direction, while they extend in- 
finitely along [011]. The reason for using Si(211) is that 
the metals are expected to nucleate at the step edges, 
and form atomic nanowires. 

The scanning tunneling microscopy (STM) images ob- 
tained by Baski et &\H for the Ga/Si(211) system showed 



the formation of atomic wires of Ga extending along the 
[011] direction with vacancies at every fifth or sixth atoms 
interval i.e., the surface structure becomes lxN with 
N=5-6. The Ga wires were separated by 9.4 A along the 
[111] direction. These authors also provided a theoretical 
model supporting their experimental observations. The 
same system was also studied by Gonzalez et alii^ and 
their STM images revealed the formation of two atomic 
chains of Ga extending along the [011] direction within 
the 9.4 A distance along the [111] direction. They also 
found that both the chains have vacancies at six atoms 
interval along the [011] direction i.e., the surface struc- 
ture is 1x6. The authors, in addition, performed a 
detailed theoretical calculations in support of their ex- 
perimental results. From the experimental and theoret- 
ical results for Ga/Si(211), it is clear that a stable un- 
interrupted Ga chain can not be formed on the Si(211) 
substrate. The In/Si(211) system has also been stud- 
ied experimentally by means of STM, low energy elec- 
tron diffraction (LEED) and Auger electron spectroscopy 
(AES). It is found that the surface pattern is 1 x 7^ They 
also presented an intuitive model in support of their re- 
sults. Wang et al>i^ performed LEED experiments on the 
Al/Si(211) at sub-monolayer coverage and they found the 
surface pattern to be 1 x 6. The question we raise here 
is whether or not a stable uninterrupted atomic Al chain 
may be formed on the Si(211) substrate. Since, no de- 
tailed electronic structure calculations have so far been 
made for the precise atomic arrangements of Al on the 
Si(211) surface, the possibility of formation of a stable 
continuous Al wire on Si(211) is also worth investigat- 
ing. In this paper we therefore perform extensive total 
energy calculations for the Al/Si(211) system. 

From theoretical perspective it is important to know 
the reconstruction of the clean surface prior to the place- 
ment of any metal on the surface. As far as the recon- 
struction of the clean Si(211) surface is concerned, there 
have been several studiesi^iSiiLiS. However, Baski et 
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FIG. 1: Atoms in perspective form the ideal Si(211) 2x4 
super-cell. The bottom layer Si atoms are passivated by hy- 
drogen atoms (small circles). Surface terrace, trench and edge 
atoms are denoted as T, Tr and E respectively. X, Y and Z 
directions correspond to [111], [Oil] and [211] respectively 

ali 9 . found (using STM) that the clean Si(211) is really 
unstable and consists of nanofacets with (111) and (337) 
orientations. It has also been shown that the Si(211) ori- 
entation is regained upon metal adsorption^i^ In view 
of this result we therefore consider the ideal bulk ter- 
minated Si(211) substrate for the Al atomic wires study. 
In our investigations, we systematically increase Al cov- 
erage on the Si(211) surface and see how the adsorption 
sites and the surface structure change as we go towards 
a monolayer (ML) coverage. We will show that atomic 
arrangement of Al on Si(211) will depend on the chem- 
ical potential of Al (/zai) which in turn depends on the 
experimental conditions such as the substrate tempera- 
ture, the Al vapor pressure in the effusion cell and the 
Al vapor pressure on the substrate. Our study reveals 
that the interrupted Al chains with period five/six (sim- 
ilar to the case of Ga) may be formed. However, most 
significant result in this investigation is that a stable un- 
interrupted Al chain extending along the [Oil] direction 
assisted by an interrupted chain can be formed at about 
1 ML coverage. 

The paper is organized as follows. In Section [n] we 
present the parameters used in the pseudopotential den- 
sity functional calculations. The results and discussions 
are presented in Section [TTT1 Finally, in section IIVI we 
summarize our principal findings. 

II. METHOD 

Total energy minimization calculations are carried 
out within the density functional theory (DFT) in con- 
junction with the pseudopotential approximation. The 
Si(211) surface is represented in a repeated slab geome- 
try. Each slab contains seven Si(211) layers with a vac- 
uum region of 12 A. For example, in the 1x4 super-cell, 
each layer contains eight Si atoms: two along [111] and 
four along [Oil]. The top layer (within the 1x4 super- 
cell) contains four edge and four terrace Si atoms. It is 
noted that in an ideal Si(211):l x 1 surface each layer 
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FIG. 2: Filled circles represent Si atoms of the bulk termi- 
nated Si(211) surface within a 1 x 4 supercell. First, second, 
and third column of atoms are terrace, trench and edge atoms 
respectively. Dangling bonds of all the atoms are shown. Var- 
ious kinds of symmetric sites are designated as G, V, F, and M 
sites respectively. The numerals 1, 2, 3, and 4 indicate identi- 
cal sites displaced along the Y direction in the super-cell. For 
example, Gl, G2, G3 and G4 (Groove sites) are identical sites 
displaced by 3.84 A along the Y direction in the super-cell. 



consists of two distinct atoms. The Si atoms in the bot- 
tom layers have their dangling bonds saturated by hydro- 
gen atoms (see Fig. ^| . Since the edge atoms have two 
dangling bonds each, the trench atoms have one dan- 
gling bond each and the terrace atoms have one dangling 
bond each, we require sixteen hydrogen atoms to satu- 
rate all the dangling bonds at the bottom of the slab in 
the 1x4 super-cell. We have used super-cells of various 
sizes (1x2, 1x3, 1x4, 1x5, 1x6, and 1x7) in our cal- 
culations. The top five Si layers are relaxed for geometry 
optimization while the two lower-most Si layers and the 
hydrogen atoms are held fixed to simulate the bulk-like 
termination. The wave functions are expanded in a plane 
wave basis set with a cutoff energy \k + G\ 2 < 250 eV. 
The Brillouin zone (BZ) integration is performed within a 
Monkhorst-Pack (MP)2£ scheme using four inequivalent 
k-points. It has been established earlier 17 that the energy 
cutoff, the set of k-points, the number of layers in the 
slab and the amount of vacuum region considered here 
give sufficiently converged values for total energies. Ionic 
potentials are represented by Vandcrbilt-type ultra-soft 
pseudopotentials^i and results are obtained using gener- 
alized gradient approximation (GGA) 22 for the exchange- 
correlation potential. Preconditioned conjugate gradient 
is used for wave function optimization and a conjugate 
gradient for ionic relaxations. The Z axis is taken per- 
pendicular to the Si(211) surface, while X and Y axis are 
along [111] and [Oil] respectively. The VASP code^ is 
used for our calculations. 
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III. RESULTS AND DISCUSSIONS 

The structures and energetics for Al on the Si(211) 
surface at 1/8, 1/4, 1/2 and beyond 1/2 ML coverages are 
systematically analyzed here to examine the formation of 
Al atomic wires on the Si(211) substrate. Note that one 
monolayer corresponds to one atom per Si atom on the 
top layer of the ideal bulk terminated surface leading to 
~ 5.56 x 10 14 atoms/cm 2 . 



A. Structure and energetics of Al at 1/8 ML 
coverage on Si(211) 

Here we discuss the structure and energetics of Al at 
1/8 ML coverage and hence we need to place one Al atom 
on the surface of 1 x 4 super-cell. At such a low coverage 
the Al atoms are practically isolated from each other on 
the surface. The Si(211) surface offers various kinds of 
sites for Al to bind. Based on the symmetry and avail- 
able dangling bonds on the surface, we consider all the 
probable binding sites which are designated as G, V, M, 
and F sites respectively in Fig. [21 We use numerals 1, 2, 
3 and 4 to label identical sites displaced by 3.84 A along 
the [Oil] direction in the super-cell. For example, Gl, 
G2, G3 and G4 (G sites) are identical sites displaced by 
3.84 A along the [Oil] direction in the super-cell. The 
binding energies (BE) of the Al atom at different kinds 
of sites are given in Table H The binding energy of the 
Al atom is defined as BE = -[E(A1 + Si) - E(Si) - E A i] 
where E(A1+Si), E(Si) and E A i are total energy of Al ad- 
sorbed super-cell, total energy of the super-cell without 
Al and atomic energy of Al respectively. 

From Table H we notice that a G site (BE=4.53 eV) 
is most favorable for Al adsorption followed by a F site 
(BE=4.45 eV), and a V site (BE=4.31 eV) respectively. 
Binding energy of Al at a V site is close to that at a M 
site. In our calculations, we let the Al atom (except at 
V sites) relax in all directions so as to reach the local 
energy minimum. For V site, the Al atom is allowed to 
relax both along Y and Z directions. It is reasonable 
that the Al atom favors to bind at a site where it can 
satisfy itself by sharing charges with three neighboring 
Si atoms. At a G site the Al atom is surrounded by two 
edge Si atoms and one trench Si atom. The distance of 
the Al atom (at a G site) from the nearest edge Si atoms 
are ~ 2.6 A and that from the nearest trench Si atom is 
also ~ 2.6 A. On the other hand, the Al atom at the F 
site also shares its charge with three Si neighbors and the 
Al-Si bonds at those sites are equally strong as that at a 
G site. For example, the distances of the Al atom at a F 
site from its neighboring terrace and two trench Si atoms 
are 2.6 A , 2.5 A and 2.5 A respectively. Furthermore, 
the charge density plot in Fig.|21clearly shows that the Al 
atom at both the G and F sites makes strong bonds with 
three neighboring Si atoms. We also notice that the edge 
Si atoms form dimers irrespective of the position of the 
Al atom at a G site or a F site. However, as the Al atom 



TABLE I: Binding energies (BE) for Al at Gl, Fl, VI and 
Ml sites on the bulk terminated Si(211) at 1/8 ML coverage. 



Site: 


Gl 


Fl 


VI 


Ml 


BE (eV): 


4.52 


4.45 


4.31 


4.30 



(a) Al at a G site (b) Al at a F site 




FIG. 3: Total charge distribution on a plane just above the 
top Si layer at 1/8 ML coverage: (a) when the Al atom is 
placed at a G site and (b) when the Al atom is placed at a 
F site. In both figures, circles in the first, second and third 
columns are terrace, trench and edge Si atoms respectively 
while the Al atoms are indicated by arrows. 



forms bonds with the neighboring Si atoms, it introduces 
strain in the underlying lattice. Therefore the net energy 
gain due to Al adsorption will also depend on the amount 
of strain introduced by the Al atom during the process 
of forming bonds with the neighboring Si atoms. From 
Fig. 3(b) it is evident that one of the edge Si dimers in 
the super-cell is distorted due to the strain introduced by 
the Al atom adsorbed at a F site. As a whole, a G site 
turns out to be most favorable. The net energy gain in 
this process is the difference between the binding energy 
of Al at a G site and its chemical potential, AE = (BE 
- \hai\) ~ 0.85 eV, where fiAl = —3.68 eV is the bulk 
chemical potential for aluminum. However, in practice, 
the energy gain will be less than 0.85 eV as the bulk 
chemical potential is the upper bound of /iAi- In other 
words, a maximum of ~ 0.85 eV is gained by adsorbing 
an Al atom at a G site. 

Our calculations show that the adsorption of Al atom 
on the terrace site is not possible because one needs to 
supply an energy of ~ 1 eV to place an Al atom on the 
terrace Si atom. We also studied the energetics for the 
replacement of terrace, trench and edge Si atoms by Al 
atoms. Our results indicate that replacement of trench 
and edge Si atoms by Al atoms are energetically not pos- 
sible, while an energy of 0.5 eV is gained due to the re- 
placement of a single terrace Si atom by an Al atom. We 
note that the chemical potential for Si (usi) is taken to 
be the bulk value ~ -5.43 eV. This information will be 
useful in the later part of our discussion at higher cover- 
ages of Al. 
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B. Structure and energetics of Al at 1/4 ML 
coverage on Si(211) 

To analyze the structure and energetics of Al on 
Si(211) at 1/4 ML coverage we place two Al atoms on the 
bulk terminated Si(211) surface of the 1x4 super-cell. 
Based on the results at 1 /8 ML coverage and on physical 
grounds, we consider various reasonable combinations of 
a couple of sites where Al atoms may prefer to bind. The 
composite sites that we consider are G1G3 (one Al atom 
is placed at Gl site and the other is placed at G3 site on 
the surface of the 1x4 super-cell), G1F3 (one Al atom 
is placed at Gl site and the other is placed at F3 site), 
G1M3 (one Al atom is placed at Gl site and the other 
is placed at M3 site), F1F3 (one atom is placed at Fl 
site and the other is placed at F3 site) and M1M3 (one 
atom is placed at Ml site and the other is placed at M3 
site) respectively. Without doing any calculations and 
just based on the results at 1/8 ML coverage, one may 
conclude that the composite site, G1G3 should be prefer- 
able for Al. However, in the presence of two Al atoms the 
surface may undergo further reconstructions due to the 
interaction with the surface Si atoms and hence, G1G3 
may not be the most favorable composite site. We there- 
fore, perform extensive calculations for Al adsorption at 
all the composite sites. The average binding energies per 
Al atom are given in Table ITU By examining Table ITTI 
it turns out that G1G3 configuration is the least favor- 
able (BE per Al is ~ 4.2 eV) configuration while the most 
favorable configuration is M1M3 (BE per Al is ~ 4.4 eV). 

We note that while M sites turned out to be least fa- 
vorable among the G, F, V and M sites at 1/8 ML cov- 
erage, it becomes most favorable at 1/4 ML coverage. 
In other words, the Al atoms move from the G sites to 
M sites through lower energy channels (with a barrier 
height of ~ 0.3 eV) as the coverage goes from 1/8 ML to 
1/4 ML. The strain induced in the underlying lattice due 
to the adsorption of Al atoms are responsible for such 
happenings. At 1/4 ML coverage, the surface consists 
of Al chains lying between terrace and trench atoms and 
extending along the [Oil] direction with an interatomic 
distance of 7.68 A. Therefore the preferable binding sites 
for Al and hence the arrangement of surface atoms are 
highly coverage dependent. 

The total charge distribution on a plane just below the 
Al layer for M1M3 and G1G3 configurations are shown 
in Fig. We observe that the Al atoms in both the 
configurations form strong bonds with their neighboring 
Si atoms. We however note that while the edge Si atoms 
in the M1M3 configuration form strong dimers (dimcr 
length 2.3 A) to reduce the total energy of the system, 
they are unable to form strong dimers (dimer length 2.9 
A) in the G1G3 configuration due to the strain induced 
by the Al atoms. This is one of the reasons for favoring 
M1M3 configuration over the G1G3 configuration. At 
this 1/4 ML Al coverage, the average binding energy per 
Al atom is slightly lower (by 0.1 eV) than that at 1/8 
ML coverage. We conclude that at 1/4 ML coverage of 



(a) M1 M3 configuration (b) G1 G3 configuration 




FIG. 4: The total charge distribution on a plane just below 
the Al layer at 1/4 ML coverage of Al on Si(211): (a) the most 
favorable configuration, M1M3 and (b) the G1G3 configura- 
tion. Al atoms are indicated by arrows. The small circles 
represent the terrace, trench and edge atoms respectively as 
one goes from left to right column. 



TABLE II: Average Binding energy (BE) per Al atom at 1/4 
ML coverage. A set of two sites combinations on the bulk 
terminated Si(211) are considered and they are denoted as 
G1G3, G1F3, G1M3, F1F3 and M1M3 respectively. 



Site: 


G1G3 


G1F3 


G1M3 


F1F3 


M1M3 


BE (eV): 


4.22 


4.26 


4.33 


4.37 


4.43 



Al the surface structure symmetry becomes 1x2 and 
the maximum energy gain per Al atom in the process of 
adsorbing two Al atoms at the Ml and M3 sites on the 
surface of the 1x4 super-cell is 0.75 eV. 



C. Structure and energetics of Al at 1/2 ML 
coverage on Si(211) 

The Al coverage is now increased to 1/2 ML i.e., four 
Al atoms are placed on the surface of the 1x4 super- 
cell. Based on previous results at low coverages and the 
dangling bonds available on the surface we consider three 
plausible combinations of four sites on the surface. They 
are designated as M1M2M3M4 (four Al atoms are placed 
at Ml, M2, M3 and M4 sites available on the surface of 
the super-cell), G1G3M2M4 (four Al atoms are placed 
at Gl, G3, M2 and M4 sites available on the surface of 
the super-cell) and G1G2G3G4 (four Al atoms are placed 
at Gl, G2, G3 and G4 sites available on the surface of 
the super-cell) respectively. The average binding energy 
per Al atom for all the configurations are given in Ta- 
ble IHII We note that average energy gain per Al atom 
for the configurations G1G3M2M4 (BE ~ 4.3 eV) and 
G1G2G3G4 (BE ~ 4.4 eV) are very close to each other. 
The charge density plot for the G1G3M2M4 configura- 
tion is shown in Fig. |Sfa). We note from Fig.FJJa) that 
all the four Al atoms make bonds with their neighboring 
Si atoms and the edge Si atoms form dimers. All the Al 
atoms remain more or less on the same plane and satu- 
rates all the surface dangling bonds. Therefore, for the 
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TABLE III: Average Binding energy (BE) per Al atom at 
1/2 ML coverage. A set of four sites combinations on the 
bulk terminated Si(211) are considered and they are denoted 
as M1M2M3M4, G1G3M2M4 and G1G2G3G4 respectively. 



Site: 


M1M2M3M4 


G1G3M2M4 


G1G2G3G4 


BE (eV): 


3.88 


4.30 


4.37 



(a) G1 G3M2M4 (b) G1 G2G3G4 




FIG. 5: The charge distribution just below the Al layer at 
1/2 ML coverage of Al on Si(211): (a) second best favorable 
configuration, G1G3M2M4 and (b) the most favorable con- 
figuration, G1G2G3G4. Al atoms are indicated by arrows in 
the figures. 



G1G3M2M4 configuration, the highly terraced Si(211) 
surface becomes more or less flat decorated by parallel 
zigzag Al chains separated by 9.4 A and extending along 
the [Oil] direction. However, for the most favorable con- 
figuration, G1G2G3G4, the terrace Si atoms remain un- 
saturated, all the Al atoms remain on the same plane 
and make strong bonds with the neighboring trench and 
edge Si atoms. Thus, at 1/2 ML Al coverage, the Si(211) 
surface becomes a flat surface decorated with straight 
parallel Al chains separated by 9.4 A and running along 
the [Oil] direction (see charge density plot in Fig.JSfb)). 
The average energy gain per Al atom in the process of 
bringing four Al atoms from the source and putting them 
at G1,G2,G3, and G4 sites is 0.7 eV which is substan- 
tial. A comparison of the configurations G1G3M2M4 and 
G1G2G3G4 suggests that the complete saturation of sur- 
face dangling bonds does not necessarily lead to the most 
favorable structure here and hence the surface strain in- 
duced by the adsorbed Al atoms becomes a significant 
factor for the most favorable structure. We therefore 
conclude that straight Al chains lying between the trench 
and edge Si atoms (i.e., on the G sites) and extending 
along the [Oil] direction can be formed. We find that 
the binding energy per Al atom decreases with increas- 
ing coverage. 




FIG. 6: The surface energies (eV) for Al chains having va- 
cancies with period four (1 x 4), five (1 x 5), six (1x6) and 
seven (1x7) relative to the no vacancy case are plotted as a 
function of /iai- Here the surface energy per lxl super-cell 
surface area (~ 36.11 x 10 _16 cm 2 ) is considered for compar- 
ison among chains with different vacancy periods. Note that 
the chemical potential of Al is given with respect to the chem- 
ical potential of bulk Al. 



D. Structure and energetics of Al beyond 1/2 ML 
coverage on Si(211) 

The formation of Al chain at 1/2 ML coverage is pos- 
sible under ideal conditions because in experiments it is 
difficult to restrict the coverage at a predetermined fixed 
value. Therefore we want to examine the evolution of 
surface structure when the Al coverage is just below 1/2 
ML and also when it is somewhat larger, tending to- 
wards 1 ML. A continuous chain lying between the trench 
and edge sites produces strain in the system which may 
be relieved by forming vacancies along the chain (like 
GaiSii&Si). We therefore compare the surface energies 
corresponding to interrupted Al chains as a function of 
the chemical potential of Al. Figure shows the varia- 
tion of surface energy as a function of hai for interrupted 
Al chains (lying on the G sites) with different vacancy 
periods. The solid line with squares correspond to the 
surface energy for the formation of an Al chain that is 
interrupted at every fourth site and resulting surface pat- 
tern becomes 1x4. Similarly the solid lines with circles 
(1x5), stars (1 x 6) and the solid line with no symbol 
(1 x 7) correspond to Al chains with vacancy periods 5,6 
and 7 respectively. The surface energies are plotted rela- 
tive to the surface energy for the continuous chain (chain 
without vacancy). We note that the surface energies for 
1x2 and 1x3 surface structures are not plotted here 
because they are very large compared to those shown in 
Fig. We know that the upper limit of //ai is the bulk 
chemical potential of Al (/xf" !fe ) which is taken as the 
reference in the figure. However, the practical value of 
Hai depends on the experimental conditions. The ther- 
modynamic relation between the Mai an( ^ ^ ne practical 
value of /iAi is given by the relation^ 

Mai = Mai* - k B Tln( Pc / Ps ) (1) 



H AI (eV) 

FIG. 7: The surface energies (eV) as a function of /j,ai for the 
situations where zero (solid line), one (solid line with inverted 
triangles), two (solid line with squares), three (solid line with 
plus signs) , four (solid line with stars) and five (solid line with 
circles) terrace Si atoms are replaced (in the 1x5 supercell) 
by Al atoms respectively in the presence of an interrupted Al 
chain with vacancy period five on the groove sites. Surface 
energies relative to that for the no vacancy case is plotted. 
Note that the chemical potential of Al is given with respect 
to the chemical potential of bulk Al. 
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FIG. 8: The surface energies (eV) as a function of /uai for 
the situations where zero (solid line), one (dotted line), two 
(solid line with diamonds), three (solid line with squares), four 
(solid line with plus signs), five (solid line with stars) and six 
(solid line with circles) terrace Si atoms are replaced (in the 
1x6 supercell) by Al atoms respectively in the presence of an 
interrupted Al chain with vacancy period six on the groove 
sites. Surface energies relative to that for the no vacancy case 
is plotted. Note that the chemical potential of Al is given 
with respect to the chemical potential of bulk Al. 



where p c is the Al vapor pressure in the effusion cell and 
p s is the Al vapor pressure at the sample and T is the 
sample temperature. Therefore, /iAi is always less than 
^Stiifc Experimental conditions can be set up to make 
this difference up to ~ 1 eV. We observe from Fig.EJthat 
the interrupted chain with a vacancy period of five is 
energetically most favorable up to 1 eV below the Mai ■ 
However, in the range 0.7 eV < |^ai| < eV, the next 
favorable chain is the one with six vacancy period. These 
interrupted Al chains with vacancy period five and six 
corresponds to an Al coverage slightly lower than 1/2 
ML. However, to find the lowest surface energy structure 
we need to go beyond 1/2 ML coverage. 

Therefore, we next consider the interrupted chain with 
vacancy period five and see how the surface energy 
changes if we add more Al atoms to the surface. Al- 
though the adsorption of Al atoms at the terrace sites is 
energetically not possible, terrace Si atoms may well be 
replaced by Al atoms. Figure [7| shows the variation of 
surface energies for replacement of zero, one, two, three, 
four and five terrace Si atoms by Al in the presence of an 
interrupted Al chain (on the groove sites with vacancy 
period five) as a function of /Uaj- The solid lines with in- 
verted triangle, square, plus sign, star, and circle symbols 
correspond to the situations where one, two, three, four 
and five terrace Si atoms are replaced by Al respectively 
in the presence of a 1 x 5 interrupted chain. The solid 
line corresponds to the no replacement case. We clearly 
notice that there are three different regimes for /zai fa- 
voring different Al configuration, (i) For, eV < |/iAi| 
0.3 eV, the replacement of all terrace Si atoms are favor- 
able in the presence of a 1 x 5 interrupted chain on the 
groove sites i.e., the Si(211) surface will consist of one in- 
terrupted Al chain with vacancy period five on the groove 



sites along with a continuous Al chain on the terrace sites 
extending along the [Oil] direction (solid line with cir- 
cles), (ii) For 0.3 eV < |^Al| ^ 0.5 eV the Si(211) surface 
favors an interrupted Al chain on the terrace sites with 
period five in the presence of the interrupted Al chain 
on the groove sites with a vacancy period five (solid line 
with stars), (in) For |/xai| ^ 0.5 the surface does favor 
any replacement of terrace Si atoms by Al atoms in the 
presence of the interrupted Al chain on the groove sites 
with a vacancy period five (solid line) i.e., the surface 
consists of only interrupted Al chains with period five on 
the groove sites. 

We already saw (in Fig. [f)J| that the interrupted Al 
chain on groove sites with the vacancy period six was the 
second most favorable structure in the practical regime 
of hai . It is therefore desirable to study the possibility of 
replacement of terrace Si atoms by Al atoms in the pres- 
ence of lx 6 interrupted chain also. Figure [S] shows the 
variation of surface energies for replacement of zero, one, 
two, three, four, five and six terrace Si atoms by Al in 
the presence of 1 x 6 interrupted Al chain (on the groove 
sites) as a function of /xai- The dotted line and the solid 
lines with diamond, square, plus, star, and circle symbols 
correspond to the situations where one, two, three, four, 
five and six terrace Si atoms are respectively replaced by 
Al atoms in the presence of 1 x 6 vacancy chain on the 
groove sites. The solid line corresponds to the no replace- 
ment case. Similar to the Fig. [7] we again find in Fig. [S] 
that there are three different regimes for /iAj favoring dif- 
ferent Al replacement/no-replacement configurations, (i) 
For, eV < |^ai| i5 0.3 eV, the replacement of all terrace 
Si atoms are favorable in the presence of 1 x 6 vacancy 
chain on the groove sites i.e., the Si(211) surface consists 
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FIG. 9: The surface energies (eV) per lxl super-cell surface 
area (~ 36.11 x 10 _16 cm 2 ) as a function of (iai for the follow- 
ing surface structures: (a) the uninterrupted Al chain on the 
terrace sites along with an interrupted chain on the groove 
sites with period five vacancy (dotted line), (b) the inter- 
rupted Al chain on terrace sites with period five vacancy along 
with an interrupted chain on the groove sites with period 
five vacancy (solid line with circles), (c) the uninterrupted 
Al chain on the terrace sites along with an interrupted chain 
on the groove sites with period six vacancy (solid line), and 
(d) the interrupted Al chain on terrace sites with period six 
vacancy along with an interrupted chain on the groove sites 
with period six vacancy (solid line with stars). Surface ener- 
gies relative to that for the no vacancy case is plotted. Note 
that the chemical potential of Al is given with respect to the 
chemical potential of bulk Al. 

of one interrupted Al chain with vacancy period six on 
the groove sites along with a continuous Al chain on the 
terrace sites extending along the [Oil] direction (solid line 
with circles), (ii) For 0.3 eV < |mai| £ 0.5 eV the Si(211) 
surface favors an interrupted Al chain on the terrace sites 
with period six in the presence of the interrupted Al chain 
on the groove sites with a vacancy period six (solid line 
with stars). However, (iii) for \/jIa\\ <1 0.5 eV the surface 
does not favor any replacement of terrace Si atoms by Al 
atoms in the presence of the interrupted Al chain on the 
groove sites with a vacancy period six(solid line) i.e., the 
surface consists of only interrupted Al chains with period 
six along the groove sites. 

Comparing Figs. EJ |H1 we can conclude that for 
I MAl | > 0.5 eV the replacement of terrace Si atoms by Al 
is energetically unfavorable. The most favorable surface 
structure is 1 x 5 where the surface consists of interrupted 
Al chains (with vacancy period five) lying between trench 
and edge Si atoms and extending along the [011] direc- 
tion. When eV < |mai| ^ 0.5 eV, Figs. □ and suggest 
that the most probable surface structure around one ML 
Al coverage can be deduced by comparing the surface en- 
ergies for four situations: (a) the uninterrupted Al chain 
on the terrace sites along with an interrupted chain on 
the groove sites with period five vacancy, (b) the inter- 
rupted Al chain on terrace sites with period five vacancy 
along with an interrupted chain on the groove sites with 
period five vacancy, (c) the uninterrupted Al chain on 



Al wires 




FIG. 10: The surface consists of two kinds of Al chains (indi- 
cated by arrows) : one of them in uninterrupted and the other 
one is interrupted with a vacancy at a period of six. Both the 
chains are extending along the Y direction. This is the most 
favorable configuration for eV < \^ai\ ^ 0.5 eV. 



the terrace sites along with an interrupted chain on the 
groove sites with period six vacancy, and (d) the inter- 
rupted Al chain on terrace sites with period six vacancy 
along with an interrupted chain on the groove sites with 
period six vacancy. These calculations present us with 
a rich set of possibilities for creating supported atomic 
wires on Si(211). Figure El shows the relative surface 
energies for such four surface structures in the range of 
|/iAi| between and 0.5 eV. For eV < \fi A i\ ^ 0.3 eV 
the surface structure corresponding to the solid line is 
most probable one (1x6) where the surface consists of 
two kinds of Al chains extending along the [011] direc- 
tion: one is continuous and the other is interrupted with 
a vacancy period six (see Fig. 1 1 1 ft . However, the sur- 
face energy represented by the dotted line is very close 
to that represented by the solid line and therefore the 
surface structure may also be 1 x 5 where the surface 
consists of two kinds of Al chains extending along the 
[111] direction: one is continuous and the other is inter- 
rupted with a vacancy period five. In any case in the 
range, eV < |/xai| iS 0.3 eV, an interrupted Al chain 
with period five or six appears followed by a continu- 
ous Al chain. The interrupted chain, relieves the strain 
from the system to enable the formation of continuous Al 
chain. In this sense, the interrupted Al chain assists the 
formation of the continuous Al chain. In other words, in 
the range, eV < \hm\ ^ 0.3 eV, continuous Al chains 
can be formed assisted by the interrupted Al chains with 
vacancy period five or six. It is worth recalling that the 
formation of continuous Ga chain is not possible on the 
Si(211) surfaced. In the range of 0.3 eV < \fi M \ < 0.5 
eV, the surface energies for the solid lines with circles 
and stars (see Fig. EJ) are lowest and they are very close 
to each other and hence we conclude that in this range 
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FIG. 11: Band structures: (a) for the surface consisting of 
continuous Al chains and interrupted Al chains with period 
1x6; (b) for the surface consisting of only interrupted Al 
chains with period 1x6. The dotted lines indicate the fermi 
levels. 

of fi ai uninterrupted Al chains will not be formed, i.e. 
the surface will consist of only interrupted chains and 
the surface structure will be either 1 x 5 or 1 x 6. Since 
the surface structure and the arrangement of Al atoms 
on the Si(211) depends on \xai which again depends on 
the experimental conditions, the observed surface struc- 
ture upon Al deposition will depend on the experimental 
conditions. In other words, a desired surface structure 
among the many (discussed above) will be obtained by 
controlling the experimental condition. For example the 
uninterrupted Al chain assisted by interrupted Al chain 
may be formed if /iai is typically maintained within the 
range of 0.3 eV from the bulk value. Our results indi- 
cate that the 1x6 surface pattern of Al on the Si(211) 
observed by Wang et al>i£ is one of the most probable 
structures. 

Band structures for two different surface structures are 
shown in Fig. 11. In Fig. 11(a), the band structure for 
the surface consisting of continuous Al chains and inter- 
rupted Al chains with period 1x6 (most favorable surface 
structure for eV < |^ai| ^0.3 eV) reveals that the sur- 
face is metallic in nature. On the other hand, Fig. 11(b), 
indicates that the band structure for the surface con- 
sisting of only interrupted Al chains with period 1x6 
(one of the favorable structures for 0.3 eV < \hai\ ^ 0.5 
eV) is semiconducting in nature. We therefore conclude 
that below 1 ML Al coverage the Si(211) surface may 
be metallic or semiconducting in nature depending on 
the experimental conditions while the Ga covered Si(211) 
surface is only semiconducting in nature^. However, like 



Ga/Si(211) an experimental study of Al/Si(211) is desir- 
able. 



IV. SUMMARY 

Energy minimization calculations are performed to 
study the surface structures, due to adsorption of Al at 
various coverages on the Si(211) surface. Well below 1 
ML coverage, we find that though the average binding 
energy per Al decreases slowly with increase of the cover- 
age, the favorable sites change dramatically. For example 
at 1/8 ML coverage the Al atoms favors to bind at the 
groove (G) sites with an average binding energy ~ 4.52 
eV and at 1/4 ML coverage Al atoms prefers to bind at M 
sites with an average binding energy ~ 4.43 eV. However, 
around 1/2 ML coverage Al atoms prefer to nucleate at 
the groove sites to form chain structure extending along 
the Y[011] direction with an interruption by a vacancy 
at regular intervals (with vacancy period 5 or 6). As 
the coverage increases further, the replacement process 
comes into play, i.e., the terrace Si atoms are replaced 
by the Al atoms to form an uninterrupted/interrupted 
Al chain along the terrace sites in addition to the pres- 
ence of the interrupted Al chain along the groove sites. A 
detailed analysis of our results reveals that for |^ai| > 0.5 
eV, only interrupted Al chain along the groove sites with 
a vacancy period of five will be formed and for 0.3 eV 

ImaiI ^ 0.5 eV interrupted Al chains with vacancy pe- 
riod five/six will be formed both along the terrace and 
groove sites. It is significant to note that for < \fi\i\ ^ 
0.3 eV, continuous Al chains along the terrace sites as- 
sisted by the interrupted Al chain (with vacancy period 
five/six) along the groove sites can be formed. We note 
that the experimentally observed 1x6 surface pattern for 
Al/Si(211) corresponds to one of the favorable structures 
presented here within the eV < \^a\\ 0.5 eV. A novel 
finding is that an uninterrupted Al chain may be formed 
on Si(211) unlike the case of Ga and In on Si(211). Fur- 
thermore, the Al covered Si(211) surface may be metallic 
or semiconducting in nature depending on the experi- 
mental conditions, while the Ga covered Si(211) surface 
is only semiconducting in nature below 1 ML coverage. 
One would expect that since Al (3s 2 3p 1 ) and Ga (4s 2 4p x ) 
belong to the same column of the periodic table, they 
should manifest similar physical and chemical properties. 
This is certainly true at a gross level but their behavior at 
surfaces can be very different due to differences in lattice 
constants but more significantly due to major differences 
in their cohesive energies. It should be noted that cohe- 
sive energy of Al (~ 4 eV) is at least ~ 0.5 eV larger (more 
strongly bound by this amount) than for Ga. When it 
comes to binding on surfaces of other materials, there is 
a competition between cohesion and bonding and hence 
it is not surprising that Al and Ga may stabilize in differ- 
ent configurations on Si(211). The STM experiments on 
the Al/Si(211) system will be helpful in revealing both 
the continuous and interrupted Al wire structures. 
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